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We study flavor violating processes in the production or decay of a neutral pseudoscalar 
meson P° in the framework of a type III two Higgs Doublet Model with four-texture Yukawa 
couplings. We use a version of the model where Yukawa interactions of neutral Higgs bosons 
allow for flavor change at the tree-level, but conserves CP symmetry. The physical Yukawa 
couplings respect CP-invariance due to the Hermiticity conditions that we impose on the 
fermion mass matrices. We focus on all possible — > / ± P° and P° — > l + l'~ decay channels, 
where I, I' arc charged leptons. We find that these processes provide complementary informa- 
tion on quark and lepton FCNC Yukawa couplings. In particular flavor violating parameters 
in the quark sector, \sb and Xdbi are significantly constrained by present experimental data, 
whereas the corresponding parameters in the leptonic sector are less constrained. 
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I. INTRODUCTION 



Flavor-changing neutral current (FCNC) processes and charge-parity (CP) symmetry violation 
observed in the quark sector are known to be duly explained in the standard model (SM) framework 
via the quark mixing mechanism [1]. Nowadays, experimental searches of CP violation and FCNC's 
in the leptonic sector are among the most interesting problems in particle physics since they would 
eventually shed some light on the origin of mixing and masses of leptons. Actually, since LF 
violating processes originating from mixing of neutrinos are expected to occur at an unobservable 
small level 

as 

, their observation would clearly indicate evidence for New Physics beyond the SM 
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In this paper, we study the LF violating processes that may occur in the production or decay of 
a neutral pseudoscalar meson P° due to the exchange of a neutral pseudoscalar Higgs boson (A ). 
To be more specific, we will focus on LF violation in all possible — > l^P° and P° — ► decay 
channels (I, I' are charged leptons) that are allowed by kinematics. Although our main focus is on 
LF violating processes, for the purposes of comparison we will also consider LF conserving decays 
of neutral heavy mesons by assuming that they are induced by the Higgs boson A alone. 

With the advent of B meson factories, large sets of r lepton pairs have been accumulated by 
both BaBar and Belle experiments Qj allowing to improve previous bounds on flavor violating 
r and B meson decays. We note also that the flavor violating decays of our interest have been 
studied in a large variety of New Physics models and using different approximations [f], Q]. A 
comprehensive set of bounds on R-parity violating couplings from these LFV decays were obtained 
in 8|. In this paper we study these decays in the framework of a general two Higgs doublet model 
type-Ill (2HDM-III) [9( which contains flavor violation via Yukawa interactions of the neutral 
Higgs bosons at the tree-level. In this paper we work in a type III 2HDM with four-texture 
Yukawa couplings, where the Yukawa matrices are assumed to be Hermitean, this implies that CP 
symmetry is respected by the Yukawa interactions of neutral Higgs bosons which largely simplifies 
the analysis since only the pseudoscalar Higgs boson contributes to the decay amplitudes of the 
processes under consideration. It is worth noticing that a restricted set of these r decays were 
considered in Ref.[lo| in the 2HDM-III but without assuming CP conservation. 

This paper is organized as follows. In section 2, we provide the Yukawa interactions for the 
neutral Higgs bosons. In section 3, we derive the LFV effective four-fermion hamiltonian and we 
provide the numerical values of couplings relevant for our calculations. The formulae needed to 
evaluate the two-body decay rates are given in section 4; from a comparison with the experimental 
upper limits on branching fractions we derive the bounds on flavor violating couplings. Finally, in 
section 5, we present our conclusions. 



II. FLAVOR VIOLATION IN THE 2HDM-III 

All lepton flavor violating processes in the SM have vanishing decay probabilities at the tree 
level and even at 1-loop level they are completely negligible [2j and beyond the reach of present 
and planned experiments. Therefore, flavor violation at a level accessible to experiments can be 
expected only in the framework of New Physics scenarios. In this section we discuss the Yukawa 
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couplings of fermions in the 2HDM-III which can induce flavor-changing neutral current processes 
both in the leptonic and quark sectors. We also provide numerical values of couplings that enter 
the hadronic matrix elements. 

In the next subsections we shall derive the Yukawa interactions of the scalar (hP, H°) and 
pseudoscalar (A ) neutral Higgs bosons with the c harg ed lepton and quark sectors of the 2HDM- 



III by assuming four-texture mass matrices 
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12, 



ia|. 



A. The charged lepton sector 



The Yukawa lagrangian of the 2HDM-III for the charged lepton sector is given by: 



C Y = yLl&iIrj + YLLi^lRj + h.c 



(1) 



where = (</>i~2> ^i 2) T denote the Higgs doublets and Li denote the doublet of left-handed 
leptons. The specific choices for the Yukawa matrices Y\ 2 define the versions of the 2HDM known 
as type I, II or III. 

After spontaneous symmetry breaking the charged lepton mass matrix is given by, 
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Mi = -f=(v 1 Yi + v 2 Yj), 



(2) 



Now we assume that both Yukawa matrices Yj and Y% have the four-texture form and are Hermitian; 



following the conventions of 



1 



the lepton mass matrix is then written as: 



Ci 
C? Bi B l 



(3) 



B* A t ) 



such that when B\ — > one recovers the six-texture form. We will also consider the following 
hierarchy, | Ai | 3> | B[ | , | Bi | , | Ci | , which is supported by the observed fermion masses in the 
SM. 

Because of the Hermiticity condition, both B[ and A\ are real parameters, while the phases 
&B,c °f Ci and B\, can be removed from the mass matrix Mi by defining: Mi = P^MiP, where 
P = diag[l, e'* c , e^* s+ * c - ) ], and the mass matrix Mi includes only the real parts of Mi. The 
diagonalization of M is then obtained by means of an orthogonal matrix O, such that the diagonal 
mass matrix is Mi = O t MiO. 
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The lagrangian (pQ) can be expanded in terms of the mass-eigenstates for the neutral (hP, H°, A ) 
and charged Higgs bosons (H^). The interactions of the neutral Higgs bosons are then given by, 



L. Y 



g ( raj 
2 \mw 



U 



cos a v2 sin(a 
-da + 



0) ( m w 



cos j3 



g cos /3 



Y 



2'J 



+ 



+ i 



sjnci . \/2 cos(o! — /3) /mjy 

mi 



cos /? 

— tan (35ij + 



5 cos (3 
y2 Z' 



mi 



2*7 



2*j 



7 5 A° ) lj + h.c. 



(4) 



g cos /3 

where a denotes the mixing angle which is used to define the physical mass eigenstates hP and 
H in the CP-even Higgs sector, and tan j3 is the ratio of vacuum expectation values of the two 
Higgs doublets [lj|. In the above expression, the terms proportional to 5ij corresponds to the 
modification of the 2HDM-II over the SM result, while the terms proportional to ft denotes the 
new contribution from the 2HDM-III. Thus, the physical fermion-Higgs couplings respect CP- 
invariance, despite the fact that the Yukawa matrices include complex phases; this follows because 
of the Hermiticity conditions imposed on both Y[ and Y\ matrices. 

The correction terms to the lepton flavor conserving (LFC) and violating (LFV) couplings, 
depend on the rotated matrix ft = T PYkP^O. We can derive ft, by assuming that Y 2 has a 



four-texture form, namely: 



yl 



( 



\ 



M | > I B, 



2 > 



B 



2 h 



Co 



(5) 



C 2 
C| B2 B2 
B* 2 A 2 ) 

Since the orthogonal matrix that diagonalizes the real matrix M\ with the four-texture form, has 
the form: 



O 



( I A 2 A 3 (A-A!) / AiA 3 (A 2 -A) 

' V A (X 2 ~X 1 ){Xj^XTj A{\ 2 -\-,){\ z -\ 2 ) 
J X^-^Aj— 
V \ (A 2 -Ai)(A 3 -Ai) 
. / A 1 (A-A 2 )(A-A^" 
\ A(A2-Ai)(A 3 -Ai) 



A 2 (A-A 2 ) 
(A2-A1XA3-A2) 
A 2 (A-A!)(A 3 -A) 
A(A 2 -Ai)(A 3 -A 2 ) 



AiA 2 (A-A 3 ) \ 
A(A 3 -A!)(A 3 ^A27 
A 3 (A 3 -A) 



(A 3 -Ai)(A 3 -A 2 ) 
A 3 (A-A!)(A-A 2 ) 
A(A 3 -Ai)(A 3 -A 2 ) 



(6) 



where m e = mi =| Ai \,mu = 777,2 ■ 
of ft acquires the general form, 

Yl 



m 3 



A3 I and T) = X2/1TI2, then the rotated form 



T PY 2 l P ] 








' (ft)n (ft 


)12 


(ft 




(Y 2 l hi (ft 


)22 


(ft 


)23 


{ (fthi (ft 


).32 


(ft 


>33 / 



(7) 
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The full expressions for the resulting entries of this matrix have a complicated form. To derive 
a convenient approximation, we will consider that the elements of the Yukawa matrix Y 2 l exhibit 
the same hierarchy as the full mass matrix, namely: 



C 2 = C2x l^pL , ( 8 ) 



B 2 = b 2 ^/(A- A 2 )(m 3 -A) , (9) 
B 2 = b 2 (m 3 -A + X 2 ) (10) 
A 2 = a 2 A. (11) 

In order to keep the same hierarchy for the elements of the mass matrix, it is found that A must 
satisfy the condition (7713 — m 2 ) < A < m 3 . Thus, we propose the following relation for A: 

A = m 3 (l-f3z), (12) 

where z = m 2 /m 3 <C 1 and < (3 < 1. 

If we now introduce the dimensionless matrix x, we get: 



>ij V 



Xij e 11 . (13) 



161 ] not only because of the appearance 



This expression differs from the usual Cheng-Sher Ansatz 
of the complex phases, but also in the form of the real parts Xij = \Xij\- 
If we expand in powers of z, the entries of the matrix \ become: 

Xil = [h -(4e i ^+c 2 e-^)]7 7 + [a 2 + b 2 -(b* 2 e i ^+b 2 e- i ^))P 

X12 = (c 2 e~^ -b^-^ + h-ib*^ +b 2 e-^)](3] 

X13 = {a 2 -b 2 eT^ B )r}\ffi 

X22 = b 2V + [a 2 + b 2 - {b* 2 ^ B + b 2 e-^ B )]P 

X23 = (b 2 e~^ B -a 2 )^ 

X33 = a 2 (14) 

As it can be easily checked, the diagonal elements xu are rea l an d the phases appear in the 
off-diagonal elements; we note that these elements are not constrained by current low-energy phe- 
nomena. Furthermore the LFV couplings satisfy some relations, such as: | X23 1 = | X13 1 ? which 
simplifies somehow the number of free parameters. 
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Finally, in order to perform phenomenological studies it is found convenient to rewrite the 
lagrangian given in Eq. ^ in terms of the Xi/s as follows: 



-y 



to,- 



cos a sin(a - 0) ( Jm^mj 

6 v + —S 7T — J Xij 



+ 9 -k 
2 

+ 2 1 



mw ) cos 
rrij \ sina 



y/2 cos/3 



mw 



mw J cos 

mi 



cos(a - 0) ( Jmirrij 
da H ;= — I | Xij 



m w 



tan 5ij + 



y/2 cos0 
1 

y/2 cos0 



mw 
jmirrij 

m w 



Ijh 



Xij 



j 5 LA + h.c, 



where, unlike the Cheng-Sher Ansatz, Xij (* / j) are now complex parameters. 



B. The quark sector 



(15) 



The Yukawa lagrangian for the quark fields in the 2HDM-III model is written as: 

C q Y = Y?Q L $ 1UR + Y 2 U Q L <S> 2 UR + YfQL^dR + Y 2 d Q L $ 2 d R + h.c, (16) 

where Ql denote the doublet of left-handed quarks. Again, the specific choices for the Yukawa 
matrices Y q 2 (q = u, d) define the versions of the 2HDM known as type I, II or III. 

Following a similar procedure as in the case of the charged leptons we can derive the Yukawa 
couplings for the physical neutral Higgs bosons (hP,H°,A°). If we introduce the dimensionless 
matrix Xqq 1 '- 



Y q 

1 2 



/ 



m q m q i 



qq' 



Xqq' 



qTTlqi il& 
-Xqq'd 



(17) 



to replace the Yukawa matrices in the rotated basis, we get the explicit form of the Yukawa couplings 
for the neutral Higgs bosons: 



L,y 
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a ° dd ' H 7T ' Xdd ' 

J cos y/2 cos 



mw 



+ - d 



, % 9 A 
+ J d 

+ 2 U 

+ 2 U 

, { 9 - 
~2 U 



mw 

m d \ sma cos(a — 0) ( Jm d m d i 
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y/2 cos/3 
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^/m d m dl 

mw 



Xdd 



, _ sin(a - 0) f yjm u m u , ' 

mw J sin yJ2 sin \ mw 

m u \ cos a cos(a — 0) ( yjm u m u i 

Ouu' 7^ — ; — ( — | Xuu' 



mw 



J sin/3 



y/2 sin/3 



m — ) cot 5 uu i + 
mw 



y/2 sin/3 



mw 
y/m u m u , 

mw 



X'b 



d'H° 
d'h° 
"f 5 d'A° 
u'H° 



u'h° 



j 5 u'A + h.c, 



(18) 
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where d (respectively, u) stands for down-type (up-type) quarks. 

In our discussion above, the flavor changing couplings are contained in the entries Xff> these 
dimensionless parameters are expected to be of order one or below and we will consider that the 
bounds obtained from data are significant if this turns out to be the case. In our calculations below, 
we will neglect the contributions of x// in the flavor conserving vertices. Note that processes with 
flavor violation at both vertices, like r — > fiKs or B° — > e/i, are allowed to occur in our model. 
However, we will not consider this type of decays because data allows to bound the product of 
flavor violating couplings and not these couplings themselves. 

It is interesting to note that in the above version of the 2HDM-III, the assumption that the 
Yukawa matrices are Hermitean leads to CP conservation, and this implies that the scalar (pseu- 
doscalar) Higgs bosons couple only to scalar (pseudoscalar) fermionic currents. This will greatly 
simplify the analysis of this paper since only the pseudoscalar Higgs boson will contribute to the 
decays of our interest. 



III. EFFECTIVE HAMILTONIAN AND COUPLING CONSTANTS 

The Feynman graphs corresponding to the two-body decays of our interest are mediated by a 
pseudoscalar Higgs boson ^4° and are shown in Figure 1. 





(b) 



FIG. 1: Two-body decays mediated by the pseudoscalar Higgs boson A : (a) tau lepton decays and, (b) 
Pseudoscalar meson decays. 



In the local four-fermion approximation, which is valid for the energy scales involved in the 
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decays under consideration, it is more convenient to work with the effective hamiltonian given by 



H = V2G1 



w 



rn 



A" 



1,1' 



+ h.c. 



(19) 



where the coupling constants gAff can be easily identified from the interaction lagrangians given 
in eqs. (1) and (2). The above interaction hamiltonian allows us to cover all the cases for flavor 
violation at the lepton and quark vertices as studied this paper. 

The relevant hadronic matrix elements that will enter in our calculations are the following: 

(ny + m q ) (0\iq'j 5 q\P) = h% , (20) 

where q, q' denote the quarks (q' 7^ q only when flavor violation occurs at the meson vertex). Note 
that the h q p constants can be directly related to the usual pseudoscalar meson decay constants fp 
(namely, hp = fpmp for all the cases except the 7r , 77, rf mesons). 

In this work we will use the isospin limit where m u = = 3.3 MeV. For the other quark 
masses we will use [r?| : m s = 0.104, m c = 1.27, mj, = 4.2 GeV. In our numerical evaluations, we 
will use the following values of the decay constants of the unflavored mesons [la ] (here q = u,d): 

2.38 x 10~ 3 GeV 3 , 
2.0 x 10~ 3 GeV 3 , 
1.6 x 10~ 3 GeV 3 , 
-53 x 10~ 3 GeV 3 , 



hi 



K 

h q 



h s - 
""q 

h s v , = 65 x 10~ 3 GeV 3 . 

Equivalently, we can use information on the pseudoscalar decay constants from Refs. 

fa = (130 ± 5) MeV, 

f K0 = (155.5 ± 0.9) MeV, 

f D o = (205.8 ± 8.9) MeV, 

f B o = (216 ± 22) MeV, 

f Ba = (260 ± 29) MeV . 

The quark content of unflavored mesons are the following: 

7T° = —= (uu — dd) 
v2 

77 = 77s cos 9p — 771 sin#p 
77 = 7/8 sin#p + 7/1 cos 9p 



(21) 



13, 



19 



20 



2li: 



(22) 



(23) 
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where 778,1 denote the octet and singlet isoscalar mesons and 9p —20° is the mixing angle of 
pseudoscalar mesons. 

IV. FLAVOR VIOLATION IN PRODUCTION/DECAY OF P° MESONS 

In this section we study the coupling of a pseudoscalar meson P° to a leptonic neutral current. 
We consider that the flavor change can occur at either the hadronic or leptonic vertex. When we 
compare our results with the current experimental bounds on branching fractions, we derive the 
constraints on the relevant flavor changing Yukawa couplings of the A Higgs boson. 

The decay widths in the case of r lepton decays are given by: 



r 2 



8tt \m A 



2 2 , \ l l 2 {m 2 T ,m 2 ,m 2 p ) 
mi) — 777.pl — 



X9Ari 



P\^2gAq Q >q'i5q\o 



q,q 



while the corresponding widths for meson decays are: 

i2 / \ 4 



x 9au> 



(0\J29A qq 'q'l 5 q\P) 



q,q 



In the above expressions we have defined \(x, y, z) = x + y + z 2 — 2xy — 2xz — 2yz. 



(24) 



(25) 



The hadronic matrix elements required for our numerical evaluations are the following: 

-^-c(gAuu - 9 Add) 



(0\O\n) = 


i 

V2 (m 
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i 
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(0\O\B s ) = 


-igAsb 



(9Auu + 9 Add) 



{9Auu + 9 Add) 



hi 



{m u + m d ) 
(m u + ?77 d ) 



+ V2~9 



Ass 



2m, 



h s , 

+ V2g As " 



' 2m s 



fDm 



D 



777 c + m u 



fBmj 



m b + m d 
2 

B s 



fB s m] 



(26) 

777 b + m s 

where the operator O is defined by O = ^ 9Aijqjl§qi and the values of quark masses and decay 
constants were given in the previous section. 
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§Aij is the coupling of the Higgs neutral boson to a fermion pair ij and in our case are given as 
follows: 



9AW 



2 



mi 



tan [3 5w + 



y/2 cos/3 



my/ 



Xiv 



(27) 



and 



9Ai 
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9 Add' 
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y/2 sin/3 
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where I, I' = e, fj,, r; u, u' = u, c, i; and d,d' = d, s, b. 



^Jm u m ul 
mw 



X% 



y/m d m dl 
\/2 cos (3 V m w 



Xdd' 



(28) 



(29) 



In order to derive the bounds on the New Physics couplings, we use the upper limits on the 
branching ratios reported by the PDG [l^]. In Tables [U [Hi and llTl| we show our results as upper 
bounds on the product of leptonic and quark couplings normalized to the square of the Higgs 
boson A for r lepton and pseudoscalar meson decays. These expressions look rather long, but 
they can be useful to easily implement future updates. Tables HI and ITT1 contain the upper bounds 
on LFV couplings, while Table [nil refers to the bounds on the flavor changing quark couplings. A 
comparison of results from Tables I and II shows that experimental data on r lepton decays highly 
constrain the tI LF couplings. Conversely we observe that decays of light mesons provide very 
poor constraints on [ie LFV Yukawa couplings, despite the tightly upper limits on their branching 
ratios . 

Using the definitions of the effective couplings gAff an d the approximations noticed in section 
II we can translate the bounds shown in Tables I— III on the dimensionless couplings Xff- As 
we mentioned before, we will say that a constraint is significant if these dimensionless couplings 
turn out to be smaller that unity. In Figures [21 [3l H] we show the bounds on the Yukawa couplings 
Xer = \Xer\, X^r = \x^r\, Xefi = {Xefil, as a function of tan by taking m A o = 300 GeV, respectively. 
As we have mentioned in section II, our imposition of the Hermiticity conditions on the fermion 
mass matrices which implies CP-conservation in the Yukawa sector allows us to derive bounds 
directly on the couplings of the pseudoscalar Higgs boson. As we observe in Figs. [2HU these 
bounds become more restrictive when tan (3 increases. Conversely, these bounds are relaxed for 
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Process 


BR (PDG 2008) [17] 


n 

r — > 7r e 


< 8.0 x 10~ 8 


T — ► 7T^|X 




t~ — » ?7e _ 


< 9.2 x 10~ 8 


r - -> ?7/i~ 


< 6.5 x 10~ 8 


r~ — > f/e~ 


< 1.6 x 10~ 7 




< 1.3 x 10~ 7 



Upper limit 
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V2 (m u +m d ) UMuu 
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V2 




1 


m A0 


V2 
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m io 
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1 



( m 


u+m d ) 
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(m 


u+m d ) 




h". 




n 




u+m d ) 

h", 

n 



- gAdd) 
~ gAdd) 
+ V2gAs 

+ ^9Ass^ : 

h a , 
v 

1 2m B 



' 2»i, 



VZgAs 



-V2g As8 ^- 



< 2.14 x 10" 8 

< 2.67 x 10~ 8 

< 4.25 x 10~ 8 

< 2.27 x 10~ 8 

< 4.24 x 10~ 8 

< 4.19 x 10~ 8 



TABLE I: LFV tau decays to a light pseudoscalar mesons. 



Process 


BR (PDG 2008)[17] Upper limit 


tt° -> n+er 
tt° -» n~e + 
rj — > /i ± e =F 
rj' — > /r^e^ 


< 3.8 x 10- 10 

m A 

< 3.4 x 10~ 9 
< 6 .0xl0- 6 
<47xl0 - 4 


T^T^rfc)^-^) <4-55xl0- 2 

m^ V2 K+md) ^ A "« < 0.14 

[Wuu + gAdd) (m J„ ld) + V2.9a s , 2m " s ] < 14.81 

Ann + gAdd)j^^ + V2 9Ass ^ < 1.21 X 10 3 



TABLE II: LFV decays of light pseudoscalar mesons. 



increasing values of the ^4° Higgs boson mass because these bounds are proportional to m 2 A0 . It is 
shown in Fig. [2] that the most restrictive bound on Xer is obtained from t~ — > jfe~ decay, where 
we obtain Xer < 134 for m^o = 300 GeV and tan/3 = 50. From Fig. [3] we observe that the most 
restrictive bound on Xfir is obtained from r~ — ► j]\x~ decay, with Xfir < 6.68 for mjo = 300 GeV 
and tan/3 = 50. We recall here that Xe 



Xut in the model under consideration [ill. \\\. 



Process 



BR(PDG 2008) [17] 



Upper limit 



B° 
B° 
B° 
B° s 
B° s 
D° 
D° 



T + T - 



< 1.3 x 10- 7 

< 1.5 x 10~ 8 

< 4.1 x 10~ 3 

< 5,4 x 10~ 5 

< 4.7 x 10~ 8 

< 1.2 x 10~ 6 

< 1.3 x 10~ 6 



I 9Af±t± 



gAdb 



T^gAdb 
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9Att 
I SAee 
I 9Af±n 



m 



■gAdb 

gAsb 
gAsb 

~gAuc 



9A,i,i „ 



m b +m d 
fBm 2 B 
mb+md 
fE>m 2 B 
m b +m d 
fB s rn 2 Bs 
m b +m s 
fB 3 in 2 Bs 

m c +m u 
m c +m u 



< 9.17 x 10" 10 

< 3.12 x lO" 10 
| < 2.20 x 10~ 7 

| < 1.86 x 10~ 8 

< 5.49 x lO" 10 
| < 2.56 x 10~ 8 

I < 2.68 x 10~ 8 



TABLE III: Flavor violation in decays of pseudoscalar mesons. 
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FIG. 3: Upper bound on Ximt as a function of tan/3 for m A a — 300 GeV (See Table HJ). 



In Figures [SHU we show the bounds on the Yukawa couplings Xsb = \Xsb\-, Xdb = \Xdb\i 
Xuc = \Xuc\ as a function of tan/3 by taking m A a = 300 GeV, respectively. As it can be seen in 
Fig. [5J the most restrictive bound on \sb is obtained from B® — > decay, where we obtain 

Xsb < 1-40 x 10~ 2 for m A o = 300 GeV and tan/3 = 50. We observe in Fig. O that the most 
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10 12 




tan(3 



FIG. 4: Upper bound on Xeu as a function of tan/3 for m^a = 300 GeV (See Table ITT|) . 

restrictive bound on Xdb is obtained from B° — > decay, we obtain Xd& < 4.01 x 10~ 2 for 

m^o = 300 GeV and tan (3 = 50. We recall here that \db = Xsb in the model under consideration. 
Using the same input parameters for tan /3 and m^o , we observe in Fig. [7] that the most restrictive 
bound on Xuc is obtained from D° — > / u _ / u + decay, namely x«c < O(10 3 ) which is indeed very 
poor. We should note that the upper bounds on the Xff quark couplings are rather conservative 
as long as the SM will also give a contribution via the usual quark mixing mechanism. 



121 ] we have reported bounds on 



In the context of the same model used in this article, in Ref. 
the Xsb coupling, which are of the same order of magnitude than those presented here. However, 
we want to point out that the upper bounds reported in the present work have been obtained by 
using processes which involve only ^4°-exchange contributions, while those reported in Ref. 12] 



have been gotten by using processes which involve /i°-exchange contributions. Hence, our bounds 



121 ] depend on both parameters, tan/3 



depend only on tan/3, while the bounds reported in Ref. 
and a. 

Finally, we can compare processes involving flavor violation at one and the two vertices by 
taking the ratio of the R a , = T(P° -> /+/ / ")/r(P° -> decay rates (here I / I'). This ratio 

is independent of the hadronic parameters and of the ^4° boson mass, actually Ru> ~ (m;//2m/) • 
(x//'/sm/3) 2 for heavy meson decays into light leptons. Thus, for large enough values of tan/3 
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(typically tan/3 > 5) and under the assumption that xiV — lj R/j,e gets suppressed by a least a 
factor of the e//i mass ratio. 

V. CONCLUSIONS 

In this paper we have studied the lepton flavor violation induced by the Yukawa couplings 
of neutral pseudoscalar Higgs boson ^4° of the 2HDM-III in the two-body decays of r leptons 
and pseudoscalar mesons. Under the assumption that the Yukawa matrices are Hermitean which 
implies that CP is conserved by these interactions, we are able to get bounds on the flavor-violating 
Yukawa couplings of the pseudoscalar Higgs boson. Using present data we have found the strongest 
bounds for the lepton-flavor violating couplings in the case of r lepton decays involving (77, rf) 
mesons. Clearly, improved experimental upper limits on these decays by one or two orders of 
magnitude will produce significant bounds on the LF violating couplings. We have also considered 
LF conserving (but FCNC quark coupling) decays of neutral mesons. In this case we are able 
to find very significant constraints on the flavor-changing quark couplings under the conservative 
assumption that the SM contribution is negligible. In particular, we find Xsb < 1-40 x 10~ 2 and 
Xdb < 4.01 x 10~ 2 for the typical values m^o = 300 GeV and tan/3 = 50. We thus conclude that 
neutral pseudoscalar mesons produced in r lepton two-body decays and their decays to charged 
lepton pairs, provide information on the flavor-violating couplings that is complementary to the 
ones from the CP-even Higgs bosons. 
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